The spin glass transition of the three dimensional Heisenberg spin glass. 
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It is shown, by means of Monte Carlo simulation and Finite Size Scaling analysis, that the 
Heisenberg spin glass undergoes a finite-temperature phase transition in three dimensions. There is 
a single critical temperature, at which both a spin glass and a chiral glass orderings develop. The 
Monte Carlo algorithm, adapted from lattice gauge theory simulations, makes possible to thermalize 
lattices of size L — 32, larger than in any previous spin glass simulation in three dimensions. High 
accuracy is reached thanks to the use of the Marenostrum supercomputer. The large range of system 
sizes studied allow us to consider scaling corrections. 

PACS numbers: 75.50.Lk 75.40.Mg. 64.60.Fr, 05.50.+q, 



Recently, decisive evidence for a spin glass transi- 
tion Q in three dimensional Ising spin glasses was found 
by us |2j and by Palassini and Caracciolo [3j . These works 
demonstrated the applicability to spin-glasses of our ap- 
proach 0, [j| to Finite Size Scaling (FSS) at the critical 
temperature, as well as that of Caracciolo and cowork- 
ers 6] for the paramagnetic state (see also 0,0)- 

However, the situation is still confusing for Heisenberg 
spin glasses [1 El El E3, Q EE IIS E3 , which is the 
more experimentally interesting case (see e.g. 0, E3) 
Indeed, numerical studies are the only theoretical tool 
to achieve progress in three dimensions. Early simula- 
tions |IJ could not reach low enough temperatures due 
to the dramatic dynamical arrest of the algorithms avail- 
able at the time, and concluded that the critical temper- 
ature, T c , was strictly zero. Yet, recent studies at lower 
temperatures |16L Il7| found indirect indications of a spin 
glass transition for Heisenberg spin glasses. Matters are 
further complicated by the Villain et al. ,12j suggestion 
of a low temperature chiral glass phase, with an Ising 
like ordering due to the handedness of the non-collincar 
spin structures (see definitions below). The simulations 
of Kawamura and coworkers |l3l gave ample support 
to this spin-chirality decoupling scenario (i.e. T c = for 
the spin glass, but T c > for the chiral glass ordering). 

In order to clarify the situation for Heisenberg and XY 
spin glasses in D = 3, Young and Lee 01 have recently 
tried our FSS methods at T C J3|3|. Although parallel 
tempering only allowed them [15j to thermalize systems 
of size up to L — 12, very clear results were reached 
for the XY spin glasses. The finite-lattice correlation 
length [2fl |. was analyzed for several system sizes 
L. As expected jH El 0, the dimensionless ratio £l/L 
crosses neatly at the same T c , for the chiral glass and 
the spin glass ordering, for XY spin glasses. In the more 
important case of Heisenberg spin glasses, their results, 
although inconclusive, were interpreted also as lack of 
spin-chirality decoupling. This conclusion has been crit- 
icized by Kawamura and Hukushima |l4j |. that studied 



somehow larger systems on very few samples. 

Here we show that a finite-temperature spin glass tran- 
sition occurs for the Heisenberg spin glass in D = 3. The 
critical temperature for the spin glass transition coincides 
with that of the chiral glass. Our results rely on Monte 
Carlo simulation and FSS analysis at T r la U, |5|. We 
adapt a lattice gauge theory algorithm |2l| to our prob- 
lem. Our algorithm thermalizes L = 32 systems, well 
beyond any previous spin glass simulation in D = 3. The 
use of Marenostrum, one of the World largest computing 
facilities, during (the equivalent of) 2.6 x 10 4 Pentium 
IV computing days allowed us to simulate 4000 samples. 
The phase transition seems to be of Kosterlitz-Thouless 
type |22j . although we may not exclude a lower critical 
dimension barely smaller than three. 

We consider the Edwards-Anderson model. The 
Heisenberg spins Si = (Sf,Sj,Sf), Si - Si = 1, live on 
the nodes a cubic lattice of size L, with periodic bound- 
ary conditions. Spins interact via the Hamiltonian ((ij) 
indicates sum over all pairs of lattice nearest neighbors): 



H 



(i,3) 



Si • S*j 



(1) 



The Jii&re Gaussian distributed quenched random cou- 
plings Their mean value is zero, while their variance 
sets the energy unit. For any quantity, O, we first cal- 
culate the thermal average for the given couplings, (0}j. 
The average over the Jij, (O)j, is only taken afterward. 

In order to detect non-planar spin structures, one 
forms the chiral (pseudovector) density field [l3j |: 

C v = S l+Sti ■ (§i x §i_ &ft ) (fx = x,y,z), (2) 

where e M is the unit lattice vector along the /i axis. 

We introduce real replicas Q: pairs of spin configura- 
tions, Sj^ and S^ independently evolving with the same 
set of couplings, {Jij}, and at the same temperature, T. 
We have as well a replicated field for the chiral densi- 
ties, {21 . From the real replicas we form the tensor field 
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g^vi^i) and the chiral vector field Ki tll : (/i, v = x, y, z) 



Their Fourier transforms, 
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yield the spin glass (SG) wave vector dependent suscep- 
tibility and the chiral glass (CG) one: 
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Both for the SG and the CG case, we have [2C 
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(2%/L, 0, 0) or permutations. One has £cg 
or ^|l G for e M • fc m i n = or 2ir/L, respectively Rota- 
tional invariance at T c implies = £cg ^ or l ar S e ^- 

Model was simulated with a mixture of heat bath 
and (microcanonical) overrelaxation taken from lattice 
QC D |21l but also effective for frustrated spin mod- 
els |23|. We straightforwardly modified the implemen- 
tation for Jij = 1 in The Elementary Monte Carlo 
Step (EMCS) consists on a sequential heat bath sweep, 
followed by (the integer part of) 5L/4 sequential over- 
relaxation sweeps, to let the microcanonical wave run 
over the system. Overrelaxation fastly evolves chirality 
(in D = 1 it inverts the local chirality). One heat bath 
update is roughly as CPU time consuming as 7 overre- 
laxations. The merits of the mixed algorithm can be as- 
sessed from Fig. ^ (standard Parallel Tempering reached 
only L = 12 for the same model QJI). Overrelaxation 
may be combined with Parallel Tempering, but this is 
unnecessary at T c . Lattices L = 4, 6, 8, 12, 16, 24 and 
32 were simulated close to T c « 0.16(2) [ll [l3| (see 
table P). We extrapolate to nearby temperatures using 
bias-corrected @ data reweighting |24j . 

Thermalization is a major issue in SG simulations. Our 
thermalization tests included the by now standard log 2 
data binning (i.e. average over all samples the second 
half of the generated data, and compare this average 
with that of the second fourth of the Monte Carlo his- 
tory, the average over the second eighth, and so on) , find- 
ing compatibility for the last three bins. Another strong 
thermalization test is the consistency of the reweighting 
extrapolation j24|. As Fig. [3] shows, the reweighting ex- 
trapolation is satisfactory for our data. It is a nice check, 
because simulations at different T are completely inde- 
pendent. 

We use the quotients FSS method 4] . The used scaling 
variable is £l/L rather than the unknown L/£oo or 
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FIG. 1: Comparison of heat bath vs. heat bath plus overre- 
laxation evolution of xsg(0) ( xsg (Aw™) in the inset) from a 
hot start (EMCS: either 1 heat bath step followed by 20 over- 
relaxation sweeps or 4 heat bath sweeps). Points in the plot 
are the average over 2500 samples of 15000 successive EMCS. 
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TABLE I: Details of simulation. For each lattice size, we 
give the studied temperatures, the number of simulated sam- 
ples (in the same order) , and the number of elementary Monte 
Carlo steps per sample. In each sample we took 10 measure- 
ments, excepting for L = 32 (20000 measurements). 



(T-T C )L 1 / V . For an observable O, diverging in the large 
L limit as (O) j L=OQ oc \T — T c \ x ° , we compare (O)j in 
lattices L\ and L 2 , at the crossing temperature T il,L2 
such that £ Ll (T Lx ' h2 )/Li = £ L2 (T L i> L *)/L 2 : 



{0 Ll (T L ^))j 



L2 
Li 



(7) 



where v is the correlation length critical exponent (the 
dots stand for scaling corrections 0). The advantages 
of Eq.0 are many 0,0, 13. S3- It is easy to use. Arbi- 
trary choices on the temperature range are avoided. The 
statistical error estimation is crystal clear. One directly 
observes scaling corrections by increasing L\ and L 2 - 

Should T Ll ' L2 be obtained from £l,cg or from £l,sg? 
There is a consensus on the di verg ence of £cg at T c , while 
that of £sg is controversial [LJ, |lj|. Moreover, £l,cg 
suffers smaller scaling corrections (Fig. 0J. So, the re- 
sults in table HT1 were obtained with £l.cg- As for £l.sg 5 



3 



we treat it here as an observable scaling as . Some 
groups HI [3 E3 expect = 1, while others [H Q 
believe that = 0. 

Several features are salient in tabled (i) Once we ne- 
glect L = 4, the L evolution of our estimates is monotonic 
with increasing L x and L 2 (we agree with Refs. 0, 
within their L window), (ii) The 1/v from (9t£sg and 
<9t£cg is compatible for L > 12 . (iii) When L\ and L2 
increases, \jv systematically decreases (in a Kosterlitz- 
Thouless scenario it tends to zero) . (iv) Exponent 7sg Iv 
stabilizes for L > 12 in a typical D = 3 value Q, hence: 



^ e 1.93±0.02 

XSG OC £ CG 



(8) 



The scenario of Refs. 0,0], where £cg diverges at T c , 
while £sg (and then xsg) does not, becomes untenable. 

At criticality (table [n|, £sg scales as L 9 rather than 
L , immediately suggesting the presence of logarithmic 
corrections to scaling. Without analytical guidance, the 
detailed numerical study of logarithmic scaling correc- 
tions [23 is out of reach. Thus, we borrow the leading 
scaling corrections from the D = 2 XY model [i(J (/^ is a 
smooth scaling function, while A is a scaling amplitude): 



£l,SG 



m^cc/L) 1+ 
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logi 
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Indeed, see Fig. |21 Eq.© accounts for the L 9 scaling 
(A = 1.3, roughly 4 times that of the 2D-XY model [H). 
One also expects [2^ multiplicative logarithms in Eq. © • 
Note however that Eq.© works in a limited range of L 
for any system barely above its lower critical dimension 
since it is the ui — > limit of the standard leading cor- 
rection to scaling, L" " w (e.g. the 2 + e expansion |z7lj 
yields w « e for the Heisenberg ferromagnet, whose lower 
critical dimension is two). 

We now study £sg /L as a function of temperature (Fig. 
|3J). If we do not correct for scaling corrections (Fig. 
top), indeed all crossings are in the rather large range 
T c = 0.16(2) [H[3. However, we observe a net shift 
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of the crossings to lower T when the system sizes grow, 
(for L = 24 and L — 32 a crossing is not found on the 
simulated T range). Yet, Fig.EJ-bottom, dividing out the 
logarithm scaling corrections in Eq. © , the curves for the 
largest systems merge around T = 0.146, as expected for 
a Kosterlitz-Thouless transition (see e.g. Ref. 0). 

For the sake of completeness, let us consider our data 
for £l,sg under the spin-chirality decoupling assumption. 
In this scenario, £oo,sg diver ges at T = 0, as £00, sg = 
A/T v * {v* = 2.2, according to[l3). FSS predicts that 
£l,sg/L is a smooth function of £/£oo,sg- Take now from 
Fig. 0-top the temperature where £,l,sg/L reaches (say) 
0.44 for L = 24 and L = 32, namely T 24 w 0.1508 and 
T32 « 0.1474. In other words, £00, SG increases a factor 
32/24 w 1.33 while the temperature merely decreases 
by a 2%. Matching this with an algebraic divergence at 
T = 0, the nonsensical result v* rj 12 is obtained. 

Consider now £ CG /Z, shown in Fig. Although less 
than for £sg, the crossings shift to lower T for larger L. 
For L = 8, 12 and 16 we resolve a crossing at T w 0.155. 
For L > 16 the curves merge at T w 0.147, suggest- 
ing again a Kosterlitz-Thouless scaling, but with smaller 
scaling corrections. The L = 32 data at low T are com- 
patible with, but above, the L = 24 data. We suspect 
of a statistical fluctuation (the three L = 32 simulations 
are independent), but a crossing may not be discarded. 

In summary, we have shown that the Heisenberg 
spin glass undergoes a spin glass transition in D = 3, 
by means of Monte Carlo simulation and FSS analy- 
sis 0, 0, 0- We have adapted a lattice gauge theory 
algorithm [2l|, that thermalizes L = 32 systems, well 
beyond any previous spin glass simulation in D = 3. 
Furthermore, a large number of samples were studied in 
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TABLE II: Critical exponents from Eq.Q (T Ll,L2 from 
ei jGG ). The O were: Csg {</>■ £sg ~ L*), Ot^l/L and 
X (j/v). The Li = 24, L2 = 32 crossing is unclear (Fig. Hj, 
so we quote results at T Ll ' i2 = 0.146 (A) and 0.145 (B). 



FIG. 2: (color online) Top: For our largest systems £sg scales 
as L ' 9 , using £cg/L as scaling variable. Bottom: The re- 
sults in the upper part may be interpreted as the logarithmic 
corrections to scaling in Eq. particularly for L > 12. 
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FIG. 3: (color online) Top: £sg in units of L, versus T, for 
several lattice sizes. Lines join data obtained from a single 
simulation [24| (all simulations independent). For the sake of 
clarity we only represent data for L > 8. Bottom: as top, 
correcting the suspected logarithmic corrections, Eq.©. 
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FIG. 4: (color online) Parallel chiral correlation length in 
units of L, as a function of T, for several lattice sizes. 



Marenostrum. For studies in the low temperature phase, 
we suggest to combine our algorithm with Parallel Tem- 
pering. The large range of system sizes studied allows us 
to conclude that, at criticality, the spin glass suscepti- 
bility scales with the chiral correlation length, Eq.(JHJ). 
Therefore, a single phase transition is present in this 
problem: spin-chirality decoupling is ruled out. We ob- 
serve logarithmic corrections to scaling , explaining why 
previous investigations [1 O El [3 Q EE ES E3 
were inconclusive. Our results are compatible with 
Kosterlitz-Thouless scaling, but they are typical as well 
of a system with a lower critical dimension barely smaller 
than three. 
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